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Abstract 

The results of the entropy generation and availability energy loss analysis under conditions of oscillating 
pressure and oscillating helium gas flow in two Massachusetts Institute of Technology (MIT) test rigs - piston- 
cylinder and piston-cylinder-heat exchanger - are presented. Two solution domains, the gas spring (“single-space”) 
in the piston-cylinder test rig and the gas spring + heat exchanger (“two-space”) in the piston-cylinder-heat 
exchanger test rig are of interest. Sage and CFD-ACE+ commercial numerical codes are used to obtain 1-D and 2-D 
computer models, respectively, of each of the two solution domains and to simulate the oscillating gas flow and heat 
transfer effects in these domains. Second law analysis is used to characterize the entropy generation and availability 
energy losses inside the two solution domains. Internal and external entropy generation and availability energy loss 
results predicted by Sage and CFD-ACE+ are compared. Thermodynamic loss analysis of simple systems such as 
the MIT test rigs are often useful to understand some important features of complex pattern forming processes in 
more complex systems like the Stirling engine. This study is aimed at improving numerical codes for the prediction 
of thermodynamic losses via the development of a loss post-processor. The incorporation of loss post-processors in 
Stirling engine numerical codes will facilitate Stirling engine performance optimization. Loss analysis using 
entropy-generation rates due to heat and fluid flow is a relatively new technique for assessing component 
performance. It offers a deep insight into the flow phenomena, allows a more exact calculation of losses than is 
possible with traditional means involving the application of loss correlations and provides an effective tool for 
improving component and overall system performance. 


Nomenclature 


Latin 

letters 

A 

Control surface 

AE 

Available energy 

AEfric 

AE(due to friction) 

AEQw 

AE(due to surface heat flow) 

AEQ X 

AE(axial heat flow) 

AEDiscr Discrepancy (AE int vs. AE ex t) 

C p 

Specific heat at const, pressure 

C v 

Specific heat at const, volume 

D 

Cylinder bore (diameter) (mm) 

E 

Total energy 

h 

Enthalpy 

k 

Thermal conductivity 

ke 

Kinetic energy 

M 

Total mass; Number of inlets 

lit 

Mass flow rate 

N 

Number of outlets 

n 

Unit normal vector 

P 

Pressure 

pe 

Potential energy 

Q 

Heat transfer 

q 

Heat flux vector 

R 

Gas constant 

r v 

Volume ratio 

S 

Piston stroke; Total enthalpy 

Q 

‘-'gen 

Entropy generation 

s'" 

gen 

Volumetric entropy generation 


s Mass specific entropy 

T Temperature 

#tspc Number of time steps per cycle 
t Time (sec.) 

A t Time step (sec.) 

u Internal energy 

V Control volume 

v Mass specific volume 

W Work 

x,y,z Cartesian coordinates 
Xi Mass fraction 

Greek letters 

V Gradient 

A Change 

T| Efficiency 

jl Dynamic viscosity 

O Dissipation function 

p Fluid density 

0 Crank angle 

Subscripts 

cond Conductive 

e Exit 

gen Generation 


H High 

i Inlet; Index; Irreversible 
int Internal 

L Low 

mix Mixing 

o Outflow 

p Pressure 

r Region, Reversible 

rev Reversible 

th Thermal 

v Volume 

vise Viscous 

Abbreviations 

BDC Bottom Dead Center 
CFD Computational Fluid Dynamics 
DTF Database Template File 
GRC Glenn Research Center 
MIT Mass. Institute of Technology 
NASA Natl. Aero.& Space Admin. 
n-D (n = 1,2,3) Dimension 
RMS Root Mean Square 
SRG Stirling Radioisotope Generator 
TDC Top Dead Center 
TFE Teflon 
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I. 


Introduction 


A. Background 

Under the auspices of NASA's nuclear systems program, the Department of Energy, Lockheed Martin 
Corporation of Valley Forge, PA, INFINIA Corp. (formerly Stirling Technology Company, Inc.) of Kennewick, 
WA and NASA Glenn Research Center are developing a high-efficiency Stirling Radioisotope Generator (SRG) for 
possible use on future NASA Space Science Missions. The SRG is being developed for multi-mission use (e.g., in 
environments with and without atmospheres), including providing electric power for unmanned Mars rovers and 
deep space missions. Performance and mass improvement goals have been established for second- and third- 
generation Stirling radioisotope power systems. Multiple efforts are underway to achieve these goals, both in-house 
at GRC and under various grants and contracts. These efforts include development of a multi-dimensional Stirling 
computational fluid dynamics code, including second law analysis post-processing to separate various 
thermodynamic losses; high-temperature materials; advanced controllers; an end-to-end system dynamics model; 
low- vibration techniques; advanced regenerators; and a lightweight converter 33 . 

This paper reports the results of a second law analysis post-processing effort to evaluate various thermodynamic 
losses in the flow spaces in two Massachusetts Institute of Technology (MIT) test rigs - piston-cylinder and piston- 
cylinder-heat exchanger - operating under conditions of oscillating pressure and oscillating flow. Extensive 
numerical analyses of the surface heat transfer, pressure, flow and temperature variations and hysteresis loss in the 
flow spaces have been undertaken in our earlier study 13 and the results were shown to exhibit good qualitative 
similarities with some disagreement in magnitudes when compared with relevant experimental data and numerical 
results from the literature 20,25,31 . These numerical results are post processed for thermodynamic losses in this study. 

Heat transfer effects in the presence of oscillating pressure and oscillating flow conditions are known to be 
important but poorly understood in a wide array of technically important mechanical power systems, such as gas 
compressors, automotive engines and Stirling engines to mention a few. A major design objective for these 
reciprocating engines is to achieve good heat transfer to and from the working fluid. Until recently the heat transfer 
estimates used in these equipment have been based on extrapolation of expressions that were developed for steady- 
pressure, steady-flow conditions. Because the ordinary, steady-state convective model contains no term to account 
for the oscillatory effect in variable volumes, designs of these equipment are nowhere near optimum. Using the two 
MIT test rigs Kornhauser 20 confirmed experimentally that there was a fundamental difference between steady and 
oscillatory flow heat transfer in the variable cylinder volumes of the MIT test rigs. Building on the work done by 
Smith, Lee, Faulkner, and Chafe 8,15,21 , Kornhauser started the development of expressions suitable for oscillating 
pressure and oscillating flow conditions. He closely integrated experiment with analysis to achieve useful results 
with good qualitative but limited quantitative success. 

During the heat to useful work conversion process in reciprocating heat engines, various thermodynamic loss 
mechanisms (pressure drop, finite heat transfer, friction, mixing, leakages, etc.), due to the non-ideal nature of the 
device, influence engine thermodynamic performance. Thus another design objective for reciprocating engines is to 
keep thermodynamic losses low. This objective however conflicts with the heat transfer objective because, for 
example, heat exchangers with higher heat transfer effectiveness generally have higher flow losses 30 . It is therefore 
valuable for the designer to know the trade-off between thermodynamic losses and heat transfer. 

The second law of thermodynamics has proven to be a useful tool in identifying the mechanisms and system 
components that are responsible for thermodynamic losses and for indicating how to minimize these losses in 
practical equipment in order to improve performance. Five main approaches, developed over the last several decades 
and reported in the literature, form the basis of second-law analysis. These are, in order of decreasing 
comprehensiveness, exergy, physical-exergy, exergy-consumption, negentropy and entropy analyses. These second 
law analysis types can be grouped into two broad categories. The first (exergy and physical-exergy analyses) is more 
detailed, and examines systems in an analogous way to energy analysis. The second (exergy-consumption, 
negentropy and entropy analyses) is concerned mainly with losses due to internal irreversibilities. Exergy analysis is 
a thermodynamic approach used for analyzing and improving the efficiency of chemical and thermal processes. In 
the analysis and design of energy systems, techniques are often used which combine scientific disciplines (mainly 
thermodynamics) with economic disciplines (mainly cost accounting) to achieve optimum designs. Many 
researchers 14,26,35,36 have developed methods of performing economic analysis based on exergy, which are referred to 
by such names as thermo -economics, second-law costing and exergo-economics. The potential usefulness of exergy 
in addressing and solving environment problems have also been demonstrated by several researchers 6,7,9,19,27,38 . An 
understanding of the relations between exergy and the environment may reveal the underlying fundamental patterns 
and forces affecting changes in the environment, and help researchers to deal better with environmental damage. 
Physical exergy includes all components of exergy except chemical exergy. Thus physical-exergy analysis may be 
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sufficient if chemical exergy is not important to the analysis or if the process involves only fixed composition flow 
(e.g. steam power cycles). Physical-exergy cannot be rationally applied to processes in which compositions of 
streams change (e.g. combustion). Exergy consumption is the difference in exergy between input and output plus 
accumulations (if any) of general work potential for a system. Entropy generation destroys the available work of a 
system and is associated with irreversibilities related to three fluid dynamic processes 11 : heat transfer across a non- 
zero temperature difference, viscous friction leading to pressure drop and non-equilibrium processes such as mixing 
of different streams of fluid, unrestrained expansion, non-quasi equilibrium compression and expansion, chemical 
reaction, leakages and other nonidealities within systems. The main difference between entropy and negentropy 
analyses 18 is that negentropy is regarded as a valuable commodity while entropy is not. Except for entropy analysis, 
all other main second-law analyses adopt the negentropy viewpoint. The entropy analysis concept is employed in 
this study. 

A valuable tool in the design of a high performance engine is a numerical code which can closely simulate the 
heat transfer, fluid dynamic and thermodynamic processes inside the engine and can account for all the 
thermodynamic losses. Sage and CFD-ACE+ commercial numerical codes are employed for this numerical study. 
The Sage code is a 1-D, multi-variable thermodynamic modeling package that supports simulation and optimization 
of spring-mass-damper resonant systems and Stirling-cycle coolers and engines. The CFD-ACE+ is an unstructured 
polyhedral cell flow solver for 2-D and 3-D problems. It provides an integrated grid generation module (CFD- 
GEOM), a graphical user interface for preparation of the model (CFD-GUI), a computational solver for performing 
the simulation (CFD-ACE), and an interactive visualization program for examination and analysis of the simulation 
results (CFD-VIEW). It is integrated with a wide variety of other physics modules making it the core of a multi- 
dimensional, multi-disciplinary analysis environment. 

B. Scope/Objective 

This research effort can be broadly classified as: 

• Modeling of the MIT Test Rigs 

• Second Law Analysis via Entropy Generation 

• Second Law Analysis Post Processing of Numerical Simulation Results 

The Sage and CFD-ACE+ commercial numerical codes are employed to obtain 1-D and 2-D computer models, 
respectively, of each of the two solution domains of interest in Kornhauser’s MIT test rigs - the gas spring (“single- 
space”) in the piston-cylinder test rig and the gas spring + heat exchanger (“two-space”) in the piston-cylinder-heat 
exchanger test rig. Hereafter, gas spring and “single-space” will be used interchangeably; so also will gas spring + 
heat exchanger and “two-space”. 

Second law analyses focusing on entropy generation and its minimization have been playing a dominant role in 
recent times to understand the irreversibility in applied engineering and transport processes 2,3,12 . Bejan 1 and 
Gedeon 16 have focused on the different mechanisms behind entropy generation in applied thermal engineering. 
Numerous other investigations (mostly numerical) 23,24,28 have been performed to determine the entropy generation 
and irreversibility profiles for different geometric configurations, flow situations, and thermal boundary conditions. 

Sage and CFD-ACE+ simulation results of temperature, velocity, pressure and surface heat transfer were post- 
processed for major thermodynamic losses using entropy generation analysis. An important and primary objective of 
this study is the development and application of a thermodynamic loss post-processor to characterize the major 
thermodynamic losses inside the two MIT test rigs. The experience gained from thermodynamic loss analysis of 
simple systems such as the MIT test rigs are often useful for application to more complex systems like the Stirling 
engine. Thus this study is motivated by the future need to incorporate a thermodynamic loss post-processor in 
Stirling engine analysis codes in order to facilitate the improvement of engine thermodynamic performance through 
better understanding of the heat transfer and power losses 32 . 

II. Modeling of MIT Test Rigs 
A. MIT “Single Space” (Gas Spring) Test Rig 

Figure 1(a) shows a schematic representation of Kornhauser’s piston-cylinder test rig. The test section consists 
of a closed piston-cylinder device mounted on a compressor base. The gas spring is the volume of gas confined by 
the fixed cylinder and the moving piston. The adjustable cylinder head enables different volume ratios to be 
achieved. A volume ratio is defined to be the maximum cylinder volume (piston at “bottom-dead-center” (BDC)) 
over the minimum cylinder volume (piston at “top-dead-center” (TDC)). A piston stroke length of 76.2 mm (3.00 
in.) and volume ratio of 2.0 are used in this study. The compressor piston drives the piston in the 50.8 mm (2.00 in.) 
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Table 1 Gas Spring Dimensions 20 


Physical Quantity Symbol 

Value 

Cylinder Bore (Diameter) 

D 

50.80 mm (2.0 in) 

Piston Stroke 

S 

76.2 mm (3.0 in) 

Volume Ratio 

r v 

2.0 



Gas Spring (Base Line Grid Size: 19 x 13). 
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diameter cylinder. The piston is sealed to the cylinder with a buna-n O-ring located more than a stroke’s length from 
the piston face so that frictional heating of the cylinder wall would only affect results minimally. The piston top 
surface is brass, while the cylinder wall and head are made of steel. The apparatus is belt driven by a D.C. motor to 
allow speed adjustment. The apparatus fill line is a capillary tube of negligible volume. The gas spring dimensions 
are tabulated in Table 1. Model representations of Kornhauser’s closed piston-cylinder gas spring using Sage and 
CFD-ACE+ are illustrated in Figs. 1(b) and 1(c) respectively and discussed below. 

1. Sage Gas Spring Model (1-D) 

Figure 1(b) is a representation of the Sage graphical interface illustrating the 1-D model of Kornhauser’s MIT 
gas spring. Sage provides 8 menus - File, Display, Edit, Scan, Specify, Process, Options, and Help - each with self- 
explanatory commands and graphical edit windows with model-component palettes. When “show window” in the 
edit menu is activated, the highest level graphical edit window (Stirling machine window) and its model-component 
palette appear. Each model-component palette contains buttons arranged in tab-selected pages corresponding to 
different categories of model components. Figure 1(b) shows the Stirling machine window with six Stirling engine 
model component palettes - “Basic”, “Canisters”, “Heat Exchanger", “Phasr Moving Parts”, “Gt Moving Parts” and 
“Composite”. A Sage model is a collection of component parts assembled in a graphical edit window by clicking and 
dragging the components into the window and connected in a particular way to form a complete system. Model 
components are organized in a hierarchical structure. The root model component contains a number of sub- 
components which may themselves contain sub-sub-components. Child (sub) models are accessed by double- 
clicking on parent models. The labeled arrows sticking out of the sides of components are boundary connectors. 
They are joined together as matched pairs with numbers (2, 3 and 4) identifying the match. The labels are meant to 
suggest what information is transferred across the connection between components. Sage is also able to perform 
such functions as open and save files, increase/decrease connector level, copy, cut and paste components. 

The gas spring is modeled using Sage’s Stirling machine model component palettes. The variable volume gas 
space in Fig. 1(a) and the charge pressure inside it are modeled in the Sage code using a generic cylinder (renamed 
“gas spring cylinder” in Fig. 1(b)) and a pressure source both obtained from the “Basic” component palette in the 
Stirling machine level. The pressure source comes with a built in steady state density connection (p st d y ) and acts as 
an isobaric gas reservoir in that the density in the gas domain adapts itself so that the time average pressure is 
continuous across the connection. The cylinder- space gas and isothermal surface, child models of the gas spring 
cylinder, are obtained from the “Gas Domain” and “Cylinder Walls” component palettes in the gas spring cylinder 
level. The arrows labeled 2 in the gas spring edit window indicate a space-time varying heat flow connection Q Gxt 
between the z-face of the cylinder- space gas and the isothermal surface. The positive-facing volume displacement 
phasor and gas charge line, both child models of the cylinder-space gas, are moved up to the Stirling machine level 
for pressure connections (P phsr ) to the constrained piston (label 3) and density connections (p st d y ) to the pressure 
source (label 4) respectively. The positive-facing volume displacement phasor represents the volume change of the 
gas space due to the motion of the piston. The piston in Fig. 1(a) is modeled using the constrained piston obtained 
from the “Phsr Moving Parts” component palette in the Stirling machine level. Its child model, the negative facing- 
area phasor, is obtained from the “Mechanical Attachment” component palette in the constrained piston level. Its 
pressure connection arrow is moved up to the Stirling machine level for connection there to the gas spring cylinder. 
Figure 1(b) shows the gas spring model components communicate with each other using density (p stdy ), heat flow 
(QGxt) and pressure (P phsr ) boundary connections. 

Once a model structure is created in Sage, numerical inputs for the model component are either specified or 
modified. Sometimes user-defined variables (special output variables) can be added to model components. Sage 
numerical simulation is initiated after initializing component and overall model parameters. Sage calculations are 
performed via its solver and optimizer routines. 

2. CFD-ACE+ Gas Spring Model (2-D). 

Figure 1(c) is a 2-D representation of the gas spring model using the CFD-ACE+ code. The CFD-ACE+ code 
exploits the symmetry of the problem domain to model only one-half of the domain as shown, thereby saving on 
computer memory space and computational time. CFD-GEOM, CFD-GUI, CFD-ACE and CFD-VIEW are the 
primary CFD-ACE+ modules used in this study. 

The CFD-GEOM module is a general purpose model builder which provides comprehensive interactive 
geometry construction and grid generation capabilities. The geometry of the gas space, the problem domain of 
interest, is represented within CFD-GEOM using fundamental geometric entities (points, lines, curves, and surfaces) 
which can be manipulated (translated, rotated, projected, split, joined, etc.) as desired. Grid generation, the process 
of discretizing the problem domain with individual cells over which the flow equations are integrated follows the 
geometric representation of the gas space. The locations of the corner points of these cells constitute the “grid” or 
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“mesh” which is stored in a data base. Three grid distribution types are available: Power Law, Hyperbolic Tangent 
and Geometric. Structured and unstructured grids and different grid sizes can be generated. A baseline structured 
grid size 19x13 is presented in Figure 1(c) for illustration. The power law grid distribution type is used to obtain a 
finer grid at the boundaries for more accurate resolution of the flow features. 

After modeling is completed, the problem is formulated using CFD-GUI, an advanced graphical user interface 
that allows complete formulation of the problem. Problem formulation involves specifying the problem type, model 
options, volume conditions, boundary conditions, initial conditions and the numerical technique needed to solve the 
flow equations. 

After problem formulation, simulation is performed using the CFD-ACE flow solver module. The DTF files 
created by the CFD-ACE flow solver are loaded directly into CFD-VIEW, an interactive graphics program with 
many tools to visualize the flow physics, animate transient data sets, as well as to extract data for post processing the 
numerical results. Tools such as point probe, line probe and strip charts are used to obtain the data of required 
parameters. 

B. MIT “Two-Space” (Gas Spring+Heat Exchanger) Test Rig 

Figure 2(a) below shows a schematic representation of Kornhauser’s piston-cylinder-heat exchanger test rig. 
The test section in this rig is essentially a modified gas spring. It consisted of a fixed cylinder head with an annular 
opening leading into an annular dead-ended heat exchanger space such that flow can continuously pass between the 
heat exchanger and cylinder as the cylinder piston expands and compresses the gas more like in a Stirling machine 
cylinder. Because of the extended heat transfer surface in the annular heat exchanger, the energy flows are more 
complicated than in a simple gas spring. 

The construction of the apparatus from the compressor base to the cylinder space just below the fixed cylinder 
head including the cylinder bore and piston stroke are consistent with those of the gas spring. The heat exchanger 
annulus is 44.5 mm (1.75 in) outside diameter and 39.4 mm (1.55 in) inside diameter, for a wall-to-wall distance of 
2.5 mm (0.10 in). The annulus is 445 mm (17.5 in) long, so that a volume ratio of 2.0 for the combined cylinder and 
heat exchanger resulted in a very small cylinder clearance volume; the cylinder-head distance is nominally 2.9 mm 
(0.11 in) at top center position. The heat exchanger entrance has the same cross-sectional dimensions as the heat 
exchanger itself, the entrance corners being as sharp as could easily be machined. The inner wall of the heat 
exchanger space is made of steel. The outer wall of the heat exchanger space is steel lined with a 1.98 mm (0.078 in) 
layer of bronze-filled TFE (Teflon). This lining was chosen for thermal properties matching those of the Pyrex glass 
substrate of the surface temperature transducers. The gas spring + heat exchanger dimensions are tabulated in Table 
2 below. Model representations of Kornhauser’s “two-space” test rig using Sage and CFD-ACE+ are illustrated in 
Figs. 2(b) and 2(c) respectively and discussed below. 

1. Sage “Two-Space” Model (1-D) 

The Sage model of the “two- space” domain is obtained by modifying the Sage gas spring model via the addition 
of the annular canister representing the heat exchanger. The generic matrix in Fig. 2(b) represents the internal 
structure of the heat exchanger and is the child model of the annular canister. It is obtained from the “Matrices” 
component palette in the annular canister level. The matrix gas and isothermal surface, child models of the generic 
matrix, are obtained from the “Gas Domain” and “Matrix Solids” component palettes in the generic matrix level. 
The porosity of the generic matrix is assumed to be 1.0, implying no porous material inside the annular canister. The 
arrows labeled 2 in the generic matrix edit form indicate a heat flow connection between the matrix gas and the 
isothermal surface. The positive gas inlet, the child model of the matrix gas, models the gas flowing from the 
annular canister into the gas spring cylinder. It is obtained from the “Charge/Inlet” component palette in the matrix 
gas level. The mass flow connection arrow is moved up to the Stirling machine level for connection there to the gas 
spring cylinder (label 6). The gas spring cylinder and constrained piston models and their child models are as 
described above for the gas spring model. The only exception is the inclusion of the negative gas inlet as a child 
model for the gas spring cylinder. It is obtained from the “Charge/Inlet” component palette in the cylinder- space gas 
level and represents the gas flowing from the gas spring cylinder into the annular canister. The mass flow connection 
arrow is moved up to the Stirling machine level for connection there to the annular canister (label 6). 

2. CFD-ACE+ “Two-Space” Model (2-D) 

As with the gas spring model, the CFD-ACE+ code exploits the symmetry of the problem domain to model only 
one-half of the domain. A model of the gas spring + heat exchanger using the CFD-ACE+ code is illustrated in Fig. 
2(c) below using a course grid (100x20). 
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Figure 2(a). Kornhauser’s MIT Cylinder + HXer Test Rig. 
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Figure 2(b). Sage Model of Kornhauser’s “Two-Space” 
Domain. 


Figure 2(c). CFD-ACE+ Model of Kornhauser’s MIT 
Two Space Domain (Grid Size: 100 x 20). 
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III. Numerical Simulation Results 


Extensive numerical simulations of the flow and heat transfer phenomena under conditions of oscillating 
pressure and oscillating fluid flow inside the “single-space” and “two-space” domains were performed using Sage 
and CFD-ACE+ codes in an earlier study 13 . Simulations using the CFD-ACE+ code were performed following 
careful optimization studies of the grid size and number of cycles of piston motion using line probe, root mean 
square (RMS) and energy conservation techniques. The Courant condition 19 ’ 34 was imposed on the choice of the 
optimum time step. The study cited above compared 1-D Sage vs. 2-D CFD-ACE+ results of energy conservation 
analyses, temperature, pressure and surface heat transfer variations, and pressure-volume diagrams. Also 2-D CFD- 
ACE+ results of surface heat flux, temperature difference between the gas temperature at the radial center of the heat 
exchanger and the heat exchanger wall, pres sure- volume diagrams, velocity vectors, temperature contours and 
hysteresis loss were compared with experimental and numerical results obtained from the literature 20,25,31 . Important 
differences and similarities were highlighted between results obtained by the 1-D Sage and the 2-D CFD-ACE+ 
codes. The 2-D CFD-ACE+ results showed good qualitative agreement with results reported in the literature for both 
the “single-space” and “two-space” test rigs. For the “two-space” test rig, there were some apparent disagreements 
in the magnitudes of the heat exchanger heat fluxes and temperature differences. For a more comprehensive 
evaluation of our numerical simulation results, the reader should consult reference [13]. Overall, the achievement of 
good energy conservation results and the close agreement with relevant experimental data and numerical results 
from the literature have enabled us to proceed to focus our current effort on second law analysis post-processing of 
our numerical results for thermodynamic losses in the MIT test rigs. 


IV. Second Law Analysis via Entropy Generation 

The computational analysis of thermo-fluid systems often makes use of the mass, momentum and energy 
transport equations. For an infinitesimal control volume these equations can be written, respectively, as 4,37 


Dp 

Dt 


- pV • v 


P 


Dv 

dt 


= Pg-Vp + V • T- 


( 1 ) 

(2) 


P 


De 

Dt 


- V q - V (pv) - w'" 


( 3 ) 


where p is density, t is time, p is pressure and T is temperature. The mass specific energy e can be decomposed into 
mass specific internal, kinetic and potential energies (e = u + V 2 V 2 + gz). The parameters q and v represent, 
respectively, the heat-flux vector and the velocity vector at points that make up the control surface. The volumetric 
work transfer rate w'" represents the contribution made by the infinitesimal control volume to the overall work 
transfer rate W delivered by the finite- size control volume. The work term accounts for the flow work, pv (v is 
specific volume), required to push the mass into or out of the control volume and the work done by viscous stresses 
(shear and dissipation). The stress vector is represented by and the gravity vector by g. 

Optimization of thermo-fluid systems requires in addition to the above conservation laws, the application of the 
second law of thermodynamics. 

For a multiport flow system’s control volume the general entropy balance which is a statement of the second 
law can be written as 4 


s 


gen 


dS 

dt 


XQ r / T , + Z n^Se - Z riljSj > 0 


( 4 ) 


The left-hand side of Eq. (4) represents the rate entropy is generated within the control volume while the right-hand 
side represents the rate of entropy accumulation inside the control volume minus entropy transfer rate into the 
control volume via heat transfer plus net entropy flow rate out of the control volume via mass flow. The second law 
stipulates that entropy generation must be non-negative in all thermo -physical processes. The parameters S, Q r , T r , 
m and s denote, respectively, the total entropy, the rate at which heat is transferred into the system, the absolute 
temperature, the mass flow rate and the mass specific entropy. The subscripts i, e and r denote, respectively, inlet 
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port, exit port and region on the system boundary across which heat is transferred. An alternate statement of Eq. (4) 
in integral form can be expressed as 4 

Sge„ = KlndV = J ^£^dV + J k q • n dA + Jpsv-ndA (5) 

V V ot A f A 

In this expression, s"' en is the volume rate of entropy generation, V the volume occupied by the open system, A is the 

control surface and n, the unit vector that is locally normal to A and points away from the system. 

The second-law statement (Eq.5) for an infinitesimal control volume is obtained by transforming the surface 
integrals in accordance with the divergence theorem and by invoking the mass-continuity equation. The resulting 
statement for the volumetric entropy generation is 


Or, 


gen 



+ 



( 6 ) 


s 


fff 

gen 



q-VT + Wq 

rj-1 2 rj-1 


(7) 


Unlike Eqs. (1-3) which require that mass, momentum and energy be conserved during a process, the second 
law (Eq.4) postulates the existence of total entropy, S, a non-conserved property of state that can be created via a 
generation or production term, S gen . Whereas the energy equation (Eq.3) representing the first law of 
thermodynamics deals with the quantity of energy and transformation of energy from one form to another, the 
second law asserts that energy varies in its quality or capacity to do useful work. Water at the top of a high cliff is a 
high quality energy source because it can be used to perform work. We can use the falling water to turn a turbine and 
produce electricity. But the high quality concentrated energy in the falling water is turned into low quality dispersed 
heat energy as it reaches the bottom of the cliff. This process of energy degradation is one way of stating the second 
law of thermodynamics. The second law is the law of energy degradation. During any chemical or physical process 
the quality or capacity of the energy to perform work is irretrievably lost. The irreversible reduction in energy 
quality described by the second law of thermodynamics fits with our common-sense notion that something is lost 
when we burn fuel. We don't lose quantity - it is conserved as the first law tells us - but the second law tells us that 
what we lose is energy quality, or the capacity to perform useful work. 

The second law analysis via entropy generation is concerned mainly with losses due to internal irreversibilities. 
As noted earlier, viscous friction, heat transfer and non-equilibrium processes are the three fundamental mechanisms 
that generate entropy within thermo-fluid systems. Increase in entropy gives a measure of the extent to which the 
energy of a system is lost or unavailable for work during a certain process. It is therefore desirable to study the 
various mechanisms responsible for entropy generation in order to minimize it in engineering equipment. Entropy 
generation may be calculated in one of two ways: by entropy flow across the external boundaries of a system 
(external entropy generation) or by entropy generated by internal processes (internal entropy generation). In 
principle, the two methods of accounting should give the same answer. Discrepancies which usually arise are often 
attributable to numerical errors (finite-difference truncation errors, round-off errors, interpolation errors, etc.). 

A. External Entropy Generation 

The cyclic integral of Eq. (5) accounts for the external entropy generation from the MIT test rig: 


^ gen, cyl (ext) f 

J^dAdt 

+ 1 

J psv • n dA dt 

Period 

A T 

Period 

A 


( 8 ) 


The integral terms are the net entropy transfer with heat and mass out of the system into the surrounding 
environment. The external entropy calculation is made just over the surface between the internal calculation domain 
and the environment. 
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B. Internal Entropy Generation 

Internal entropy generation can be accounted for by tallying up the individual entropy generations in all 
internal processes. For the MIT test rigs considered in this study, possible internal entropy generation mechanisms 
include conductive heat flow in a temperature gradient and viscous dissipation resulting in mixing loss. That is: 


S 


gen (int) 


S 


gen, cond 


+ S 


gen, vise diss 


(9) 


Consider the differential energy equation that is valid for a Newtonian fluid under very general conditions of 
unsteady, compressible, viscous, heat-conducting flow 24 


P 


Dh 

~Dt 


kV 2 T + 1^- + (D 
Dt 


( 10 ) 


Note that this equation is essentially derived from Eq. (3) by incorporating the flow work in the specific enthalpy 
term (h = u + pv) thus accounting for thermal energy being convected in and out of the control volume. The shear 
stress is eliminated in the work term (along with the kinetic and potential energies) by the introduction of the linear 
momentum equation (Eq.2) leaving only the viscous dissipation term <P . Radiation heat transfer and internal 
sources of heat that might occur during a chemical or nuclear reaction are neglected leaving only conduction heat 
transfer. The parameters k and p denote the thermal conductivity, and pressure respectively. The viscous dissipation 
function <D has the form 
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(3-D) 

( 11 ) 

(2-D) 


where p is the dynamic viscosity and u, v, and w denote Cartesian velocity components. The density p and enthalpy 
h are provided by data or algebraic expressions for state relations of thermodynamic properties e.g., p = p(p, T) and 
h = h(p, T). In principle, compressible flow calculations can be made for any fluid equation of state. But in fact most 
elementary treatments are confined to monatomic ideal gases with constant specific heats such as helium, the 
working gas in the MIT test rigs. For these gases 


P = 


P 

rt’ 


R = c p - c v = const 


( 12 ) 


The parameter R is the gas constant and c p and c v are the specific heats at constant pressure and volume respectively. 

Explicit relations for differential changes in entropy can be obtained by solving for ds in the second Tds 
equation 5 

Tds = dh - vdp (13) 

By substituting v = 1/p into Eq. (13) the differential entropy change becomes 



1 dp 

pf 


Eq. (14) can be expressed in terms of the material derivative of the entropy as 


(14) 


Ds _ p Dh 1 Dp 

P Dt _ T Dt T Dt 

If we now substitute Eq. (10) in Eq. (15) we get after simplification 

Ds kV 2 T O 

p — + — 

Dt T T 


(15) 


(16) 
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Subst. Eq. (16) in Eq. (7) gives after simplification: 
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T 


(17a) 


q • VT 
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(17b) 


The internal entropy generation due to conductive heat flow and viscous dissipation in the MIT test rig can thus be 
expressed in integral form as: 


1 gen, cyl (int) 


r r Q • VT 

f 1^2— dV dt + 


Period V 


, . o 

f J T dVdt 

Period V I 


(18) 


The first term on the right-hand side of Eq. (18) which involves the squares of the temperature gradients is the 
entropy generation due to conductive heat flow while the second term which involves the squares of the velocity 
gradients is the contribution due to viscous dissipation. The sum of these two terms seems to account for all of the 
internal thermodynamic losses. The viscous dissipation term accounts for mixing loss in the flow field. Mixing loss 
is generated by mixing of the flow field due to the shear stress between fluid layers, vortices resulting from wall 
flow separation or shed vortices in the core of the fluid. A control volume approach can be used quite successfully to 
calculate the viscous dissipation term and hence the total mixing loss. Errors in the mixing loss are likely to be 
small, since all mixing processes obey the fundamental conservation equations that also govern the behavior of the 
numerical predictions. The difficulties found in using numerical predictions to estimate the mixing losses are 
thought to be due to false entropy due to numerical dissipation 10 . Although mixing is only a relatively short-lived 
phenomenon as the flow eventually mixes out completely, it is a major contributor to the total loss. This is especially 
true for turbulent flows where the effective viscosity is very large. 

Equation (18) seems to imply that the internal thermodynamic losses can be minimized by minimizing the 
magnitude of these local temperature and velocity gradients and the sum of their squares, throughout the domain of 
interest. These measures are however difficult or impractical to implement. For example, it is not practical to reduce 
the fluid viscosity which would lead to a reduction of the velocity gradients at the wall and the boundary layer 
viscous losses. Also, it is not practical to increase the fluid thermal conductivity which would reduce the temperature 
gradients at the wall, and the conductive heat transfer loss. 

C. Availability Energy Loss 

All internal entropy generations can also be characterized in terms of availability energy loss defined as 5 : 

Availability Loss = W r - ^ 

where W r and Wj (< W r ) are the net mechanical work available from reversible and irreversible heat engines 
respectively; both obtained from the same heat input Q jn and operating between a high temperature reservoir at 
temperature T H and a low temperature reservoir at temperature T L . 


The efficiency of a heat engine is defined as 5 


Tlth = 


W, 


net, out _ | _ Qput 

Qi„ 


( 20 ) 


From this definition, the net heat outflow for the irreversible heat engine q q . should be greater than the net heat 


outflow for the reversible heat engine, Q . 

For the reversible heat engine 5 

th,rev 


^2L = 1-lb 

Qin T h 


( 21 ) 
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where the Kelvin thermodynamic temperature scale (Q in /Q out ) rev = T H /T L has been used 5 . 


Thus: 


W„ 


f 


V 



and 


% = Qir 


i_Q^ 

Q 


in J 


( 22 ) 


(23) 


Substituting Eqs. (22 and 23) in Eq. (19) gives, after simplification: 


Availability Loss 




(24) 


For steady periodic process in a closed system, ^ - o an d m = 0 . Substituting the simplified Eq. (4) in Eq. (24): 

dt 

Availability Loss = T L S gen (25) 

The availability-loss concept allows us to think about entropy generation in terms of the more concrete notion 
of lost mechanical work. A loss in availability equates to a decrease of PV power in an engine. For example, in 
turbo-machines that generate shaft power (turbines) or absorb power (pumps, compressors), the rate of power lost 
owing to irreversibilities is proportional to a loss in availability and thus to the rate of entropy generation. 

The results of the second law analysis via entropy generation can be summarized as follows: 

External Entropy Generation: 

Sge„,cyl(ex,) = j> ! ^ dAdl + j> j pSV • II dA dt (26) 

Period A A Period A 


Internal Entropy generation: 
Conductive heat flow: 


Viscous Dissipation: 


Availability Loss: 
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gen, cyl(int) (cond) 


, , q • VT , 

- f dV dt 
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Period V A 


Availability Loss 


t l s 


gen 


(27) 

(28) 
(29) 
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V. Second Law Analysis Post Processing of Numerical Simulation Results 

Equations (26-28) clearly show that local entropy production depends functionally on the local values of heat 
transfer rate, temperature, density, mass-specific entropy, velocity and viscous dissipation. Thus entropy generation 
can be considered a derived quantity that can be computed by post-processing experimental or numerical flow fields. 
Equation (29) shows that local entropy production is a function of availability loss and temperature. 

A. Using Sage 

Sage entropy generation results for both test rigs are calculated from Eq. (29) using availability loss results 
obtained from Sage output file viz.: 


S 


gen 


Availability Loss 


AEfric + AEQw + AEQx ± | AEDiscr| 
To 


(30) 


The parameters AEfric, AEQw, and AEQx are available energy losses due to flow friction, surface heat flow and 
axial heat flow respectively. IAEDiscrl| is the absolute value of the discrepancy between the total available energy 
loss due to internal entropy generation (AEinternal = AEfric + AEQw + AEQx) and that due to external entropy 
generation, AEexternal. The “+” is used when it is assumed that AEexternal is greater than AEinternal and the “- 
“sign is used when the reverse is the case. These assumptions are arbitrary since the relative magnitudes of 
AEexternal and AEinternal cannot be determined apriori. A discrepancy of zero implies the total available energy 
loss due to internal entropy generation is equal to that due to external entropy generation. 

B. Using CFD-ACE+ 

1. External Entropy Generation, S gen(ext) 

The “single-space” and “two-space” models each taken as a whole constitutes a closed, isothermal, non- 
adiabatic, reciprocating system. Eq. (26) is thus reduced to 

a , f q n , (31) 

S ge„(ex.) = I J— dAdt 
Period A -L 


in which case results of surface heat transfer and temperature only are post-processed. Whereas the “single-space” 
model is limited to the use of Eq. (31) because of its simple closed system construction, the “two-space” model can 
also be analyzed, by considering its separate components - the heat exchanger space and the cylinder space - as 
open systems. For open systems, Eq.(26) is applicable and CFD-ACE+ generated results of surface heat transfer, 
temperature, mass flow rate and mass specific entropy are post-processed for external entropy generation. Figure 3 
shows the “single-space” model (Fig. 3(a)) set up as a single cylinder volume condition for analysis. In order to 
simplify the analysis, the “two-space” model (Fig. 3(b)) was partitioned into four volume conditions - outer, mid and 
inner cylinder volume conditions in the cylinder space plus the heat exchanger constituting a single volume 
condition. The external entropy generation is calculated for each volume condition with the sum of the results for the 
separate volume conditions yielding the total external entropy generation for the “two- space” model. 


Cylinder 

Outer 



Heat Exchanger 


Mid 


-K 

Piston 


Piston 


Inner 


Axis of Symmetry 


Axis of Symmetry 


(a). “Single-Space” Model 


(b). “Two-Space” Model 


Figure 3. Model Set-Up for External Entropy Generation Analysis 
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Integrated surface heat transfer rates generated each time step can be stored in output data files called assembly 
files. An assembly file is enabled by first creating an input file ( “filename. fmt”), which specifies the relevant system 
bounding surfaces, in Notepad (facility for creating and editing text files using basic text formatting) and then 
choosing the “BC Integral Output” option in CFD-ACE-GUI. To calculate the external entropy generation for the 
“single-space” model using assembly file data, the integrated heat transfer rate data are exported to an Excel 
spreadsheet for calculation of the cyclic time integral value (see Eq. 31). Division of the cyclic integral value by the 
constant surface temperature completes the calculation of the external entropy generation. The same procedure 
applies for the “two-space” model if it is analyzed as a closed system. If the “two-space” model is analyzed by 
considering its separate components as open systems, as done in this study, CFD-ACE+ generated results are post- 
processed for the external entropy generation by using a user subroutine interfaced with the CFD-ACE+ solver. Note 
that post-processing of assembly file data is only possible for closed systems. User subroutines can be used to post- 
process any CFD-ACE+ generated closed or open system data. User subroutines, especially necessary for open 
systems, require an in-depth knowledge of Fortran concepts. 

1.1 User Subroutines 

User subroutines are written in FORTRAN 90 to post-process CFD-ACE+ generated surface heat transfer, 
temperature, velocity gradients, mass flow and mass specific entropy data. The discretized domain of interest 
(interior and bounding surfaces) consists of cells (interior) and short line segments called “faces” (boundary) as 
shown in Figure 4 below. 



Figure 4. Discretized computational domain showing “faces” and “cells” 

Surface heat transfer, temperature, velocity and temperature gradients, mass flow and mass specific entropy 
values can be obtained at cell or “face” centers except at the interface between volume conditions where data values 
can only be obtained at cell centers near the interface. The velocity and temperature gradient values at cell centers 
can either be obtained directly from “get gradient function” module in CFD-ACE+ or calculated using finite 
difference approximations for the non uniform grid spacing in the computational domain 29 . The summation of 
surface, volume and time values, implied by the integrals in Eqs. (26 - 28) is facilitated with the use of do loops in 
the subroutines. In order to minimize program complexity, the subroutines are written to sum up the surface and 
volume integral values at each time step only and thereafter data is exported to an Excel spreadsheet for calculation 
of the cyclic time integrals. A carefully constructed algorithm is needed for implementation of Eqs. (26-28). 
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2. Internal Entropy Generation , S aen ( int } 

The volume conditions in the model may be further refined for calculation of internal entropy generation as 
shown in Figure 5 below in order to obtain a better resolution of the losses. Equations (27) and (28) are used to 
calculate the internal entropy generation due to conductive heat flow and viscous dissipation respectively. Results of 
heat transfer, temperature and gradients of temperature and velocity are post-processed via user subroutines written 
for both MIT models. Local and global distribution of entropy generation rates due to conductive heat flow and 
viscous dissipation can be evaluated via a histogram map. 




(a). “Single-Space” Model (b). “Two-Space” Model 

Figure 5. Model Set-Up for Internal Entropy Generation Analysis 
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VI. 


Results and Discussion 


The results of the post-processing analysis are shown and discussed below for the indicated models and 
operating conditions. 

A. Gas Spring 

Figure 6 illustrates the CFD-ACE+ results of the external and internal entropy generations as functions of the 
number of cycles of piston motion in the gas spring model. A negative external entropy generation which appears to 
violate the second law is noted for the first cycle of piston motion. This amounts to an effective source term in the 
Navier-Stokes equations because of initial transients generated by the moving piston in a flow field that is initially 
not properly defined. 



Figure 6. Gas spring entropy generation vs. Cycle number 
(10 RPM, 192.4 kPa., T wall = 294 K, Grid size = 32 x 20, #tspc = 480). 

The internal entropy generation appears insensitive to the number of cycles of piston motion one cycle earlier 
than does the external entropy generation and from the third cycle onward the two methods of accounting for the 
entropy generation appear to yield essentially the same result as they should in principle. The ~ 1.3% discrepancy 
between the external and internal entropy generation results is insignificant and probably due to numerical errors for 
the internal calculations. As expected, the external entropy generation profile is the mirror image of the graph of 
j>SQ/T i n the energy conservation plot (Fig. 6(a), [13]). 

Table 3 below shows Sage and CFD-ACE+ entropy and availability energy loss results for the gas spring. All 
CFD-ACE+ results, with the exception of the external entropy generation, are less than corresponding Sage results. 
Sage results show that the primary contribution to the total availability energy loss is the energy loss to surface heat 
flow (AEQw). With no contributions from flow friction (AEfric) and axial heat flow (AEQx) and negligible 
numerical errors as reflected by a discrepancy value of 1.031E-06 (W), Sage appears to have completely satisfied 
the principle that the two methods of accounting for the entropy generation should give the same answer. CFD- 
ACE+ calculates about 1.3% difference between internal and external entropy generation values and therefore does 
not satisfy the accounting principle as well as Sage. 

The percentage differences in the internal and external availability energy loss results between Sage and CFD- 
ACE+ are on the average about a factor of 3.4 greater than the corresponding entropy generation results. The 
external results show better agreement between Sage and CFD-ACE+ than corresponding internal results by an 
average factor of about 1.7. This is not surprising, since internal entropy calculations using 1-D Sage code are 
expected to incur more errors where the flow features inside the gas spring are more two-dimensional than one- 
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Table 3. Gas Spring Entropy Generation and Availability Energy Loss (Sage vs. CFD-ACE+) 


Gas Spring 

Sage 

(T wall = 300 K) 

CFD-ACE+: V2004 (Alpha version) 

(10 RPM, 192.4 kPa., T wall = 294 K,) 

(Grid size = 32x20, #tspc = 480, Opt. Cycle = 4) 

Internal 

External 

Internal 

External 

Entropy Generation (W/K) 

• AEfric + AEQw + AEQx 

gen,int. '-p 

^ wall 

^ AEfric + AEQw + AEQx - AEDiscr 

^ gen, ext. — y 

^ wall 

s = s + s 

gen, int. gen, cond gen, vise 

S = -H— 

gen, ext J J ^ 

dt ds A 

0.001444667 

0.001444667 

0.001432516 

0.001451513 

Availability Energy Loss (W) 

AE int = AEfric+ AEQw+ AEQx 

AE ext = AEfric + AEQw + AEQx - AEDiscr 

AF — t S 

int 1 wall ° gen, int 

AF — t S 

ext 1 wall ° gen, ext 

0.433400000 

0.433400000 

0.421159704 

0.426744822 

AEfric 

0.000E+00 




AEQw 

4.334E-01 




AEQx 

0.000E+00 




AEDiscr 

1.031E-06 



1% Difference l| between Internal and External Entropy Generation (or Availability Energy Loss) 

0.0 

1.3 


1% Difference! between Sage and V2004 (Alpha version) 

Internal 

External 

Entropy Generation 

0.8 

0.5 

Availability Energy Loss 

2.9 

1.6 


dimensional such as during the expansion stroke at crank angle 0 = 10° and 180° and compression stroke at 0 
=350° and 360° when the flow fields show large recirculation zones [13]. CFD-ACE+ is more suited to handle multi- 
dimensional flow situations. 
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B. Gas Spring + Heat Exchanger 

Figures 7 and 8 below illustrate the CFD-ACE+ results of the external and internal entropy generations 
respectively for each of the “two-space” sub-domains. Figure 7(a) shows the functional dependence of entropy 
generation on the number of cycles of piston motion and Figure 7(b) shows a histogram of the entropy generation 
values characterizing the contributions of each sub-domain to the external entropy generation. Figures 8(a) and 8(b) 
illustrate corresponding CFD-ACE+ results for the internal entropy generation. 

• As with the gas spring domain, Fig. 7(a) shows the external entropy generation plots to be independent of the 
number of cycles of piston motion beyond the third cycle. Also, the total external entropy generation plot is a 
mirror image of the graph of |§q/t i n the energy conservation plot (Fig. 6(b), [13]). The negative values for the 

external entropy generation in the cylinder sub-domains should be interpreted as the cylinder acting as “entropy 
sink”. That is, these sub-domains appear to extract entropy from their surroundings (due to heat entering the 
cylinder in these sub-domains). 


External Entropy Generation from “Two-Space” Sub-domains (CFD-ACE+) 




7(b) 
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Internal Entropy Generation in “Two-Space” Sub-domains (CFD-ACE+) 
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Figures 7&8. “Two-Space” External (7a, b) and Internal (8a, b) Entropy Generation vs. Cycle Number and Sub-Domain 
(201.7 RPM, 1.008 MPa., T wall = 294 K ; Grid size = 147 x 51, #tspc = 480). 

Fig. 8(a) shows the internal entropy generation plots to be independent of the number of cycles of piston motion 
one cycle earlier than does the external entropy generation plots, as was the case in the gas spring model. It is noted 
that the heat exchanger contributes more to the total external entropy generation than it does to the total internal 
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entropy generation. The heat exchanger’s contribution exceeds the total external entropy generation by about 81% 
(see Fig. 7(b)) but falls short of the total internal entropy generation by about 14% (see Fig. 8(b)). A clear inference 
from these plots is that the major sub-domain contribution to entropy generation is from the heat exchanger. (Note: 
In Fig. 7(b) component “external” entropy generations are external to the particular component. Only the total 
external entropy generation is solely due to “external to the system” entropy generation). 

Of the three cylinder’s sub-domains, the outer cylinder’s contribution to the entropy generation is most 
significant (external or internal). When taken as a whole, the cylinder’s contribution to both external and internal 
entropy generations is minimal. 

Table 4 illustrates entropy generation and availability loss results obtained for the “two-space” domain using 
Sage and CFD-ACE+. All CFD-ACE+ results, with the exception of the external entropy generation and availability 
energy loss results for the cylinder, are greater than corresponding Sage results. Calculation of the percentage 
differences between internal and external results of entropy and availability energy loss using CFD-ACE+ are based 
on the assumption that external results are more accurate since the calculation of the integral heat transfer rates are 
more dependable. With Sage, it does not seem to matter what the basis for calculation of the percentage differences 
is. Calculation of the percentage differences between Sage and CFD-ACE+ are based on the assumption that CFD- 
ACE+ results are more accurate since CFD-ACE+ is more suited to handle multi-dimensional flow situations. 

With the exception of availability energy loss due to axial heat flow, results of losses due to flow friction and 
surface heat flow posted by Sage are greater in the heat exchanger than in the cylinder. Sage posts negligible flow 
frictional losses in the cylinder. As in the gas spring, the numerical errors reflected by the discrepancies between the 
internal and external availability energy losses indicate that Sage does a better job than CFD-ACE+ in satisfying the 
entropy generation accounting principle that the two methods of calculating for the entropy generation should give 
the same answer. This conclusion is also reflected in the results for the percentage differences between internal and 
external entropy generations (or availability energy losses) posted by Sage (~7%) and CFD-ACE+ (-12%). 

The percentage differences between Sage and CFD-ACE total internal and external entropy generations for the two- 
space test rig are 41.8 and 51.8, respectively. From Table 3, the corresponding percentage differences for the gas 
spring are 0.8 and 0.5 % respectively. The availability energy loss results show similar trend. Thus there is much 
better agreement between Sage and CFD-ACE in the entropy generation and availability energy loss calculations for 
the gas spring, than for the two-space rig. This may possibly be due to the inability of the Sage 1-D code to 
accurately account for temperature and velocity gradients caused by flow separations from walls where there are 
changes in flow area as for example, between the cylinder and the heat exchanger in the “two-space” test rig. Sage 
1-D flow is assumed to immediately adjust to the wall boundaries through all geometrical changes in area and can 
only approximately account for the effect of flow separations, vortex formation, etc. where there are changes in flow 
area. It is noted that in the gas spring, Sage and CFD-ACE+ external calculations agree more than their internal 
calculations. The reverse is noted for the “two-space” model. 

Figure 9 shows the distribution of entropy generation inside the two-space model. The impact of the entropy 
generation due to conductive heat transfer, fluid friction and mass transfer on the efficiency of the MIT models is 
quantified. Under the specified condition (201.7 RPM, 1.008 MPa., T wa n = 294), relatively higher losses are 
experienced in the heat exchanger. This is followed by losses in the outer cylinder and inner cylinder. Contributions 
from the mid-cylinder are very negligible. The heat exchanger’s conductive heat transfer has the most impact on the 
“two-space” model efficiency followed by outer and inner cylinder conductive heat transfer. Contributions to the 
entropy generation due to viscous dissipation throughout the “two-space” domain and mid-cylinder’s mass transfer 
and conductive heat transfer are very negligible. The relatively high losses in the heat exchanger support the 
extensive effort in the Stirling community to reduce the pressure drop in the regenerator which is a form of heat 
exchanger. In compact-heat-exchanger passages for example, improvements can be made in the constructive details 
of the channels (channel shape and aspect ratio, curvature of the return channels, etc.), the actual temperature 
differences between the cold and hot fluids, the surface finishing, and the type of materials used. 
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Table 4. Gas Spring + Heat Exchanger Entropy Generation and Availability Energy Loss 


Gas Spring + Heat Exchanger 

Sub-Domain 

Sage 

(T wall = 300 K) 

CFD-ACE+: V2004 (Alpha version) 

(201.7 RPM, 1.008 MPa., T wall = 294 K,) 
(Grid size = 147x46, #tspc = 480, Opt. Cycle = 6) 

Internal 

External 

Internal 

External 

Entropy Generation (KW/K) 

^ AEfric + AEQw + AEQx 


AEfric + AEQw + AEQx - AEDiscr 

C - C _i_ Q 

gen, int. gen, cond gen, vise 

^ gen, ext “ 

dtds i 

^gen.int. — rp 

wall 

8en,ext ' T wall 

Heat Exchanger 

0.00007069 

0.00006535 

0.000137 

0.000327 

Cylinder 

0.00002245 

0.00002198 

0.000023 

-0.000146 

Total 

0.00009314 

0.00008733 

0.000160 

0.000181 

Outer Cylinder 



0.0000155 

-0.000108 

Mid Cylinder 



0.0000014 

-0.000003 

Inner Cylinder 



0.0000062 

-0.000035 

Availability Energy Loss (KW) 

AE Components 

Heat 

Exchanger 

Cylinder 

Heat 

Exchanger 

Cylinder 

Heat 

Exchanger 

Cylinder 

Heat 

Exchanger 

Cylinder 

AEfric 

0.00025860 

0.00000000 















AEQw 

0.02089000 

0.00659400 















AEQx 

0.00005789 

0.00014010 















Sub-Domain 

AE int = AEfric + AEQw + AEQx 

AE ext = AEfric+AEQw+ AEQx- AEDiscr 

AF — T S 

>rvj ^ int 1 wall ° gen, int 

ath - t S 

ext wall gen, ext 

Heat Exchanger 

0.02120649 

0.01960549 

0.040295 

0.096142 

Cylinder 

0.00673410 

0.00659400 

0.006793 

-0.042866 

Total 

0.02794059 

0.02619949 

0.047088 

0.053276 

AEDiscr 

0.0017411 

0.006188 

Outer Cylinder 



0.004565 

-0.031859 

Mid Cylinder 



0.000412 

-0.000809 

Inner Cylinder 



0.001816 

-0.010198 

Sub-Domain 

1 % Difference! between Internal and External Entropy Generation (or Availability Energy Loss) 

Heat Exchanger 

8.0 

58.0 

Cylinder 

2.0 

116.0 

Total 

7.0 

12.0 

Outer Cylinder 


114.0 

Mid Cylinder 


149.0 

Inner Cylinder 


118.0 


1 % Difference! between Sage and V2004 (Alpha version) 

Internal 

External 


Entropy Generation 

AE Loss 

Entropy Generation 

AE Loss 

Heat Exchanger 

48.4 

47.4 

80.0 

79.6 

Cylinder 

2.4 

0.9 

115.1 

115.4 

Total 

41.8 

40.7 

51.8 

50.8 
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Two-Space External Entropy Distribution 
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Figures 9. External and Internal Entropy Distribution in the “Two-Space” Model 
(201.7 RPM, 1.008 MPa., T wall = 294 K ; Grid size = 147 x 51, #tspc = 480, at Opt. Cycle = 6). 

VII. Concluding Remarks 

A discussion of the computer models of the MIT gas spring and “two-space” test rigs, post-processing of the 
numerical simulation results and entropy generation/availability energy loss analysis results for the models have 
been provided. Sage and CFD-ACE+ simulation results of temperature, velocity, pressure and surface heat transfer 
were post-processed for major thermodynamic losses via entropy generation. An important and primary objective of 
this study is the development and application of a thermodynamic loss post-processor to characterize the major 
thermodynamic losses inside the two MIT test rigs. The experience gained from thermodynamic loss analysis of 
simple systems such as the MIT test rigs are often useful for application to more complex systems like the Stirling 
engine. Thus this study is motivated by the future need to incorporate a thermodynamic loss post-processor in 
Stirling engine analysis codes in order to facilitate Stirling engine performance optimization. 

The following conclusions can be drawn from the results of this study which sought to characterize the 
irreversibilities related to heat transfer, mass flow and viscous friction occurring in the MIT test rigs via entropy 
generation: 

(1) In both the “single-space” and “two-space” models the two methods of accounting for the entropy generation 
(external and internal) appear insensitive to the number of cycles of piston motion beyond the third cycle. 

(2) In both test rigs, CFD-ACE+ does not satisfy the accounting principle as well as Sage. In the “single-space” 
model the two methods yielded essentially the same result as they should in principle. Clear differences in 
results are noted for the “two- space” model. 

(3) With few exceptions, all CFD-ACE+ results of entropy generation and availability energy loss are less than 
corresponding Sage results for the “single-space” model but greater for the “two-space” model. 

(4) Sage and CFD-ACE entropy generation and availability energy loss calculations are in better agreement for the 
“single-space” model than for the “two-space” model. This may possibly be due to the fact that CFD-ACE+ is 
more suited to handle multi-dimensional flow situations. 

(5) It is noted that in the gas spring, Sage and CFD-ACE+ external calculations agree more than their internal 
calculations. The reverse is noted for the “two-space” model. 

(6) Sage results for the gas spring show that the primary contribution to the total availability energy loss is the 
energy loss to surface heat flow. 
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(7) In the “two-space” model, the heat exchanger provides the major sub-domain contribution to the entropy 
generation. Thus the heat exchanger heat transfer has the most impact on the “two-space” model efficiency 
followed by outer and inner cylinder heat transfer. Viscous dissipation throughout the entire “two- space” 
domain and mid-cylinder mass transfer and conductive heat transfer contribute minimally to the model’s 
efficiency. 

Typically neglected and often viewed as superfluous, the second law of thermodynamics remains an esoteric 
and mysterious subject 17 particularly in computational analysis of thermo-fluid systems. When properly applied, the 
second law of thermodynamics has proven to be a very powerful tool in the optimization of complex thermodynamic 
systems. Loss analysis using entropy-generation rates due to heat and fluid flow is a relatively new technique for 
assessing component performance. It offers a deep insight into the flow phenomena, allows a more exact calculation 
of losses than is possible with traditional means involving the application of loss correlations and provides an 
effective tool for improving performance. Designers will know the cumulative amount of all losses computed locally 
in the flow domain. Entropy generation maps can be produced, and designers can use them by scanning them to 
detect critical areas (locations in which entropy generation is higher than its integral average value over the entire 
flow field). By considering the local values of entropy generation rates due to thermal and viscous dissipation, 
designers can generate a thermodynamically better design by simply trying to avoid these critical areas or re- 
computing them after a design modification has been introduced to assess local and global effects of the design 
change. 

Our understanding of loss mechanisms is however far from complete. Although numerical predictions are 
valuable in predicting the heat transfer and flow structure, there are difficulties in predicting the loss accurately. This 
is due to errors in predicting the boundary layers, transition as well as due to false entropy generation due to 
numerical dissipation. This work provides a point of reference for incorporation of loss post-processors into Stirling 
engine numerical codes. The incorporation of a loss post-processor in Stirling engine numerical codes, it is believed, 
will facilitate the optimization of Stirling engine performance. 
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